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ABSTRACT 

We estimate the baryonic (stellar+cold gas) mass function of galaxies in the local Universe by assigning a 
complete sample of Two Micron All Sky Survey and Sloan Digital Sky Survey galaxies a gas fraction based on 
a statistical sample of the entire population, under the assumption of a universally-applicable stellar initial mass 
function. The baryonic mass function is similar to the stellar mass function at the high mass end, and has a 
reasonably steep faint-end slope owing to the typically high cold gas fractions and low stellar mass-to-light ratios 
characteristic of low-mass galaxies. The Schechter Function fit parameters are (jf h~ 3 =0. 0108(6) Mpc" 3 logioM , 
M* h 2 = 5.3(3) x 10 10 M Q , and a = -1.21(5), with formal error estimates given in parentheses. We show that the 
Hi and H2 mass functions derived using this indirect route are in agreement with direct estimates, validating our 
indirect method. Integrating under the baryonic mass function and incorporating all sources of uncertainty, we 
find that the baryonic (stellar+cold gas) mass density implied by this estimate is f2 C oidbaryon = 2. x 10~ 3 , or 
8^% of the Big Bang nucleosynthesis expectation. 

Subject headings: galaxies: general — galaxies: luminosity function, mass function — galaxies: stellar content 



1. INTRODUCTION 

The distribution of the mass in collapsed baryons (cold gas 
and stars) in galaxies is a fundamental prediction of galaxy for- 
mation models. Unfortunately, to date there is no robust esti- 
mate of the baryonic mass function (MF) of galaxies, leaving 
modelers with the non-trivial task of predicting stellar masses 
or, even worse, galaxy luminosities. Discrepancy between the 
model and data may indicate a problem with the predicted dis- 
tribution of galaxy baryonic masses, or could represent poorly- 
constrained star formation (SF), stellar population or dust pre- 
scriptions. In this Letter, we present a first estimate of the bary- 
onic MF of galaxies by assigning galaxies gas fractions statisti- 
cally (based on an independent sample), under the assumption 
of a universally-applicable stellar initial mass function (IMF). 1 

The time is ripe to attempt this for the first time. With the 
advent of large, relatively complete surveys, the luminosity 
function (LF) is now well-constrained in the optical and near- 
infrared (NIR) (Gardner et al. 1997; Cole et al. 2001 ; Kochanek 
et al. 2001; Norberg et al. 2002; Liske et al. 2003; Blanton et 
al. 2003; Bell et al. 2003). Furthermore, under the assump- 
tion of a universally-applicable stellar IMF, the distribution of 
stellar masses is reasonably well-constrained, with an overall 
normalization uncertainty caused by our relatively poor knowl- 
edge of the faint end slope of the IMF (Cole et al. 2001; Bell et 
al. 2003). Crucially, there are also relatively large samples of 
galaxies with /f-band data and gas masses, allowing a reason- 
ably accurate characterization of the gas mass of galaxies as a 
function of their physical parameters (Bell & de Jong 2000). 

2. THE DATA AND METHODOLOGY 

Because of the lack of a large galaxy survey with both gas 
mass and /f-band data, we take a sampling approach, analo- 
gous to that used by Loveday (2000) to estimate the /f-band 
luminosity function from a B-band limited survey. Essentially, 

'A time-varying IMF, as speculated on by Hernandez & Ferrara (2001) or 
Ferguson, Dickinson, & Papovich (2002), would invalidate this assumption. 



we estimate a stellar MF (§2. 1) and then add representative gas 
masses to each galaxy (§2.2), allowing us to estimate the distri- 
bution of galaxy baryonic masses (§3). We assume O m = 0.3, 
n A = 0.7, and H = I00h km s" 1 Mpc" 1 . 

2.1. Estimating the Stellar Mass Function 

We construct the baryonic MF using a combined sample of 
galaxies from the Two Micron All Sky Survey (2MASS; Skrut- 
skie et al. 1997) and the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). We use the SDSS early data release (EDR; 
Stoughton et al. 2002) to provide an 84% redshift complete 
r < 17.5 sample of galaxies with accurate ugriz fluxes over 
414 square degrees, which is ^10% less than the whole EDR 
imaging area because some spectroscopic plates that were not 
attempted (Stoughton et al. 2002). The 84% redshift complete- 
ness within this area is our own direct estimate based on the 
fraction of galaxies fulfilling the Strauss et al. (2002) criteria 
that have spectra, in agreement with the EDR analysis of Naka- 
mura et al. (2003). To account for light missed by the Petrosian 
magnitude estimator (Strauss et al. 2002; Blanton et al. 2003), 
we add 15% to the fluxes of galaxies morphologically classi- 
fied as early-type using the SDSS r-band concentration param- 
eter following Strateva et al. (2001). This correction produces 
only a < 5% effect on the LFs and MFs (Blanton et al. 2003; 
Bell et al. 2003). We also correct for an ~ 8% overdensity of 
galaxies in the EDR, as estimated by comparing the number of 
10 < A" < 13.5 galaxies in the EDR spectroscopic area with that 
from the sky with \b\ > 30deg. 

We use the now complete 2MASS extended source and point 
source catalogs to augment the SDSS ugriz fluxes with /f-band 
fluxes, and for extended sources /T-band half-light radii. We 
correct 2MASS /f-band fluxes to total following a comparison 
with deeper /f-band data from Loveday (2000); for extended 
sources this amounts to a 0.1 mag correction (Bell et al. 2003). 
We do not use 2MASS J or //-band data here because we can- 
not correct the magnitudes similarly. The optical and NIR mag- 
nitude zero points are accurate to ^0.05 and ^0.02 mag respec- 
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tively, and the random errors are 0.05 mag (optical) and 0.2 mag 
(NIR). 

To estimate ^-corrections, evolution corrections, and the pre- 
sent day stellar mass-to-light ratios (M/Ls), we fit the ugrizK 
observed fluxes 2 with model stellar populations. These popu- 
lations have a range of metallicities and SF histories at a given 
redshift. We use the Pegase stellar population synthesis model 
(see Fioc & Rocca-Volmerange 1997, for a description of an 
earlier version of their model) with a 'diet' Salpeter IMF (fol- 
lowing Bell & de Jong 2001) that has the same colors and lu- 
minosity as a normal Salpeter IMF, but with only 70% of the 
mass (due to a smaller number of low-mass stars). Corrections 
derived using this technique are consistent with those used by 
Blanton et al. (2003). The stellar M/Ls we derive are within 
10% of those from the spectral modeling technique of Kauff- 
mann et al. (2003), accounting for differences in IMF; the ran- 
dom and systematic uncertainties from dust and bursts of SF 
dominate, however, and are < 25% (Bell & de Jong 2001). 
This IMF is 'maximum disk', inasmuch as IMFs richer in faint 
low-mass stars over-predict the rotation velocity of Ursa Major 
Cluster galaxies with TT-band photometry and well-resolved Hi 
rotation curves. This prescription thus gives the maximum pos- 
sible stellar M/L. Naturally, a different choice of IMF allows 
lower M/Ls. For example, the popular Kennicutt or Kroupa 
IMFs have ^37% lower M/Ls than this IMF, and are thus 'sub- 
maximal' (see Bell & de Jong 2001, for more discussion of this 
point). 

We calculate LFs and stellar MFs using the V /V max formal- 
ism (Felten 1977), taking into account foreground Galactic ex- 
tinction, ^-corrections, and evolution corrections. In Bell et al. 
(2003), we match precisely published LFs; in particular, we re- 
produce the g-band and TT-band LF and luminosity densities to 
within < 10% (Blanton et al. 2003; Cole et al. 2001; Kochanek 
et al. 2001). Furthermore, this method produces accurate stel- 
lar MFs that match the estimate of Cole et al. (2001) to ~5% in 
total stellar mass density (accounting for IMF differences), but 
can do so using LFs limited by optical or NIR magnitude limits 
(Bell et al. 2003). For this Letter, we choose 11848 galaxies 
with 13<r<17.5 and g < 17.74, which ensures that we have 
accurate g — r color estimates providing a stellar M/L accuracy 
of better than 25%, while avoiding potential biases against low 
surface brightness galaxies in 2MASS (Bell et al. 2003). The 
stellar MF estimated using this technique is shown in Fig. 2, 
along with the stellar MF of Cole et al. (2001) for comparison. 
A much more extensive description of the LF and stellar MF 
construction is given by Bell et al. (2003). 

2.2. Estimating Gas Masses 

Because there are no samples of galaxies with good number 
statistics, deep optical/NIR data and gas masses, we estimate 
the gas masses of SDSS+2MASS galaxies indirectly. We use 
galaxies from Bell & de Jong (2000) with TT-band luminosi- 
ties, half-light radii and gas masses to statistically assign a gas 
mass to every SDSS+2MASS galaxy, appropriate to its TT-band 
luminosity and half-light radius. 

Fig. 1 shows the TT-band half-light radii and luminosities for 
the late-type subsample of the SDSS+2MASS galaxies (con- 
tours) and for the comparison sample of 156 galaxies with gas 
masses (filled circles) taken from Bell & de Jong (2000). We 

2 Not all galaxies have ugrizK fluxes. We have checked that missing pass- 
bands do not significantly bias the estimated ^-corrections, evolution correc- 
tions or stellar M/Ls (but do, of course, increase the random error somewhat). 
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FIG. 1. — ^-band magnitude and half-light radius for the 6573 late-type 
SDSS+2MASS galaxies in our sample (contours) and for the 156 galaxy com- 
parison sample of Bell & de Jong (2000, filled circles). 

estimate gas masses by multiplying the Hi gas mass by 1.33 to 
account for helium, and by the morphological type-dependent 
H2 to Hi ratio presented by Young & Knezek (1989). For the 
21% of the galaxies without TT-band luminosities and sizes, we 
adopt the r-band half-light radii and estimate the Zf-band lumi- 
nosity by dividing the ^-band-derived stellar mass by the model 
K-bmd stellar M/L. Our assumptions are accurate to better than 
40% in both cases, and make no difference to our results. 

Each SDSS+2MASS galaxy is assigned a gas mass from a 
galaxy in the Bell & de Jong (2000) comparison sample with 
similar half-light radius and TT-band luminosity. For galaxies 
morphologically classified as late-type we assign gas masses 
using a randomly chosen comparison galaxy that is within a 
factor of two in size, and within one magnitude in TT-band lu- 
minosity. Galaxies are morphologically classified using the r- 
band concentration parameter following Strateva et al. (2001). 
We scale the gas mass by the difference in luminosity to con- 
serve the galaxy gas fraction. As a consistency check, we also 
assign gas masses by choosing the nearest galaxy in half-light 
radius-luminosity space, and by assigning gas to galaxies of 
all morphological types (see Fig. 2). These changes make no 
appreciable difference to our results. 

3. THE BARYONIC GALAXY MASS FUNCTION 

Our baryonic galaxy MF is shown as the solid line with open 
circles in the left panel of Fig. 2. Shown for comparison is the 
stellar MF from Cole et al. (2001), and the ^-selected stellar MF 
described above and in Bell et al. (2003). The baryonic galaxy 
MF follows the stellar MF at high masses, shows a modest off- 
set at the 'knee' of the MF, and shows a reasonably steep faint- 
end slope. This behavior is expected, as low-mass field galax- 
ies tend to have high cold gas fractions and more ongoing SF 
(Bell & de Jong 2000; Kauffmann et al. 2003b), which steep- 
ens the baryonic MF compared to optical/NIR LFs that typi- 
cally have a ~ -1 .0. The Schechter Function fit parameters are 
c/)*h- 3 =0.0108(6) Mpc-MogKjM" 1 , M*h 2 = 5.3(3) x 10 10 M Q , 
and a = -1 .21(5), where the formal error estimates are in paren- 
theses. 

In the middle panel of Fig. 2, we show the effects of the 
systematic uncertainties. In particular, the dashed, dotted and 
dot-dashed lines show the effects of changing the gas mass as- 
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FIG. 2. — The baryonic mass function of galaxies. In the left hand panel, we show the stellar MF of Cole et al. (2001, naked error bars) corrected to our 'maximum 
disk' IMF, the stellar MF derived using our SDSS g-band selected sample (solid grey line), and the baryonic MF of galaxies, assuming the 'maximum disk' IMF 
(solid line with open circles and error bars, with a Schechter Fit as the thin solid line). In the middle panel, we show different versions of the baryonic MF, illustrating 
the different sources of uncertainty. The solid line with open circles again is the baryonic MF defined using the default gas mass estimation technique, the dashed 
line shows the effect of choosing the nearest galaxy on the half-light radius-luminosity plane for estimating the gas mass, the dotted line shows the effect of allowing 
SDSS+2MASS early-type galaxies to also have gas, and the two solid lines are the effect of choosing a Kennicutt or Kroupa IMF (higher line), or a M/L consistent 
with a Bottema disk (lowest line). In the right hand panel we test the gas mass estimation method by comparing the Hi and H 2 MFs predicted using this method with 
observations. The solid line with open circles and error bars denotes the match of SDSS+2MASS spirals with randomly-selected nearby galaxies, and the dotted 
and dashed lines are as in the middle panel. The solid black curve is the blind Hi MF of Rosenberg & Schneider (2002). The H2 galaxy MF is shown as the grey 
line with open circles and error bars, and the Schechter Fit to the first observational estimate of the H 2 MF (from Keres, Yun, & Young 2003) is shown as the solid 
grey line. 



signment (see the figure caption for more details); these un- 
certainties have only small effects. The two bare solid lines 
show the effect of assuming different stellar IMFs. As stated 
earlier, we adopt as our default an IMF that has the largest 
stellar M/L permitted, without over-predicting the rotation ve- 
locities of spiral galaxies in the Ursa Major Cluster (solid line 
with open circles; Bell & de Jong 2001). Yet, the stellar M/L 
may be lower than this maximal value; thus, we show two 
cases. First, we plot (upper solid line) the increasingly pop- 
ular Kroupa, Tout, & Gilmore (1993) and Kennicutt (1983) 
IMFs, which both have M/Ls of ~70% of maximal IMF (-50% 
Salpeter). The other case is an IMF that has M/Ls of only 
40% of our maximal IMF (the lowest solid line), which cor- 
responds roughly to the 63% disk velocity contribution to the 
rotation curve as argued by Bottema (1997). It is clear that, 
assuming a universal IMF, the factor-of-two uncertainty in the 
stellar M/L dominates the error budget in terms of total cold 
baryonic mass in the local Universe. Given the stellar IMF, 
the systematic uncertainties are < 25%. We estimate these un- 
certainties using different passbands to estimate stellar mass, 
using different gas mass assignment methods, and accounting 
for the effects of dust and bursts of SF on the M/Ls (Bell 
& de Jong 2001; Bell et al. 2003). The Schechter Func- 
tion fit parameters for the Kennicutt/Kroupa case and the Bot- 
tema case are: <j>* /T 3 = 0.01 16(5), 0.0142(8)Mpc" 3 logioM" 1 , 
M* h 2 = 3.78(1 1), 2.24(8) x 1O 1O M , and a = -1 .22(3), -1 .20(3). 

It is important to make sure that our statistical procedure as- 
signs gas masses consistent with the true galaxy population by 
comparing with the Hi or H2 MF of galaxies (the right hand 
panel of Fig. 2). The solid black line with open circles and error 
bars is our Hi MF derived in this way. The dashed and dotted 
lines show the effects of using the closest galaxy to estimate gas 
mass and allowing elliptical galaxies to have gas mass, respec- 
tively. For comparison, the Hi MF of the blind Hi Aricebo sur- 



vey of Rosenberg & Schneider (2002) is plotted as the smooth 
solid curve. We also show our prediction of the H2 galaxy MF 
as the lower solid grey line with open circles and error bars. For 
comparison, we plot the Schechter Fit to the first observational 
estimate of the H 2 MF (Keres, Yun, & Young 2003). Given the 
uncertainties, it is clear that our method for estimating galaxy 
gas masses reproduces the Hi and H2 galaxy MFs with aston- 
ishing precision. Therefore, the sample of galaxies that we use 
to assign gas masses is indeed reasonably representative, and 
lends considerable credibility to our estimate of the baryonic 
galaxy MF. 

4. DISCUSSION 

Even with the factor-of-two uncertainty from the contribu- 
tion of low-mass stars to the overall stellar M/L, we can still 
draw some conclusions about the local Universe. It is clear that 
the overall efficiency of galaxy formation is very low. Firstly, 
the faint end slope of the baryonic MF is — -1 .2, which is much 
shallower than the — -2 expected for the halo MF (e.g., White 
& Frenk 1991). Secondly, integrating under the MF, we de- 
rive fWdbaryon/* = 2.4^ x 10" 3 , including the IMF and 25% 
systematic stellar M/L uncertainties. Our estimate agrees well 
with the value of 2.9 ± 1.5 x 10" 3 from Fukugita et al. (1998), 
and is preferred due to our better accounting for stellar M/Ls 
compared with Fukugita et al. (1998) who use (harder to con- 
vert into stellar mass) B-band luminosity densities assuming a 
similar IMF to the maximum-disk IMF we adopt here. Tak- 
ing the value of the total baryon density from the Big Bang 
nucleosynthesis value of O'Meara et al. (2001), and assuming 
h=0.7±0.07, we find £l co id baryon /^b ~ 8^5%, where the error 
estimates account for the uncertainties in IMF, Ho, Sib, our gas 
assignment method, and the <25% uncertainties in stellar M/Ls 
from dust and bursts of SF. Our value is quite consistent with 
the low galaxy formation efficiency characteristic of most cur- 
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rent models, which have low efficiencies at the low and high- 
mass ends because of feedback from supernovae and inefficient 
gas cooling, respectively (e.g., Cole et al. 2000). 

Accounting for the possible gaseous content of elliptical 
galaxies, for sub-maximal M/Ls, and for the effects of dust 
and bursts of SF on stellar M/Ls, the universal gas fraction, 
f g = ^coid gas/^coid baryons, should lie in the range 0.2 <f g < 0.5. 
For the 'maximal' IMF, we find f g ~ 25%. Fukugita et al. 
(1998) and Keres, Yun, & Young (2003) find values of 15- 
20% when our IMF is adopted; their slightly lower determina- 
tions stem primarily from a lower estimate of Hi mass density. 
Nevertheless, all the studies agree that f g < 0.5; therefore, the 
dynamically cold baryons (i.e., the gas and stars in disks and 
spheroids) are primarily in the form of stars, even for low stel- 
lar M/Ls. 

It is well-known observationally that cluster optical/NIR LFs 
have steeper faint-end slopes than field LFs (e.g., Trentham & 
Tully 2002). Furthermore, cluster galaxies tend to have lit- 
tle ongoing SF and little gas, so that most cluster galaxies are 
star-dominated with large stellar M/Ls (e.g., Kuntschner 2000). 
Thus, the trend of increasing faint end slope with increasing 
cluster mass noted by e.g., Trentham & Tully (2002) may be 
more naturally interpreted as a constant baryonic MF, with a 
suppression of recent SF in massive clusters of galaxies. Ob- 
viously, a deeper investigation of this issue is warranted before 
speculating any further. 

5. CONCLUSIONS 

Together with the baryonic (stellar+cold gas) luminosity- 
linewidth relation (e.g., McGaugh et al. 2000; Bell & de Jong 
2001), the baryonic galaxy MF is an ideal test of models of 
galaxy formation and evolution. In this Letter, we have es- 
timated the baryonic galaxy MF in the local Universe for the 
first time assuming a universally-applicable stellar IMF. We as- 
sign gas fractions statistically to a large sample of galaxies from 
2MASS and SDSS, using a local sample with accurate TT-band 
and gas fraction data. We cross-check this statistical procedure 
against independent Hi and H2 surveys, finding excellent agree- 
ment. The baryonic MF is similar to the stellar MF at the high 
mass end (with a slightly higher density normalization), and 
has a reasonably steep faint end slope, a ~ -1.2, due to the 
typically high cold gas fractions and low stellar M/Ls of low- 
mass galaxies. Integrating under the baryonic MF, we find that 
the baryonic (stellar+cold gas) mass density implied by this es- 
timate is tt coU baryon = 2.4^} x 10" 3 , or 8^% of the Big Bang 
nucleosynthesis expectation. This clearly implies a low overall 
efficiency of galaxy formation. 
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